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Abstract 



It is shown that the van der Waals free-energy of polydisperse fluids, as 
introduced previously ( L. Bellier-Castella, H. Xu and M. Baus, J. Chem. 
Phys. 113, 8337 (2000) ), predicts that for certain thermodynamic states 
(e.g. low temperatures and large polydispersities) the ordinary two-phase 
coexistences become metastable relative to a fractionation of the system into 
three phases, reducing thereby the polydispersity of each of the coexisting 
phases. 
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I. INTRODUCTION 



A polydisperse system can be viewed as a "continuous mixture", i.e. as the limit of 
a discrete mixture whereby the number of component species tends to infinity while the 
concentration of each species tends to zero 0. Such a system serves as a good model for 
many of the complex fluids used in various industries i. Whereas many of these industrial 
fluids have badly characterized concentration distributions one can find at present, e.g. 
in the physics laboratories, several soft condensed matter systems i (e.g. colloids, liquid 
crystals, polymers, etc) with much better controled species distributions allowing hence for 
a more detailed study of the influence of the polydispersity on various physical properties 
(e.g. phase behavior, rheology, etc.). To transpose the results and techniques elaborated for 
simple fluids and their mixtures to polydisperse or continuous mixtures is however by no 
means a simple task but has nevertheless recently become an active field of research! - !. 

In the present study, a sequel to ref.8, we will be concerned with the equilibrium phase 
behavior of polydisperse fluids. To simplify the problem, we restrict ourselves to systems of 
spherical particles (e.g. colloids) exhibiting a single polydispersity, say a size-polydispersity. 
We assume moreover the initial (parent phase) size-distribution to be monomodal, i.e. cen- 
tered around a single reference species, as suitable for the polydisperse generalization of 
one-component systems. This initial size distribution will further be assumed to be fixed, 
once and for all, by the production process of, say, the colloidal particles. (Note that some 
systems, e.g. micellar solutions, exhibit a "variable" polydispersity whereby the size distribu- 
tion is allowed to adjust itself to some externally imposed conditions!.) When these colloidal 
particles interact by excluded volume repulsions and, say, van der Waals (vdW) attractions 
their thermodynamic properties can be studied on the basis of the polydisperse generaliza- 
tion of the vdW free-energy0, such as the one already used in our earlier studiesi~l. Of 
course, such a description is far from being exact, but in general the vdW-approximation 
captures, at least qualitatively, the essence of the underlying phase behavior, as can be wit- 
nessed from several previous investigations^. As a final limitation, we would like to stress 
that our study will be limited to the fluid phases of the polydisperse system, leaving aside 
whether these phases are stable or metastable with respect to the possible solid phases. 
Indeed, the inclusion of the solid phases would require a more precise specification of the 
interaction potential than is usual in a vdW- description and, above all, would add still more 
complexity to the already fairly complex problems raised by the study of phase equilibria in 
polydisperse fluids. In this way, we will be able to focus on the central difficulty resulting 
from the replacement of the algebraic equations in finite dimensional spaces, characteristic 
of the phase behavior of discrete mixtures, by the integral equations in infinite dimensional 
spaces characteristic of the continuous mixtures. 

In our previous study@, the phase behavior of the present system was already studied 
within the same vdW- approximation, but only at the level of the two-phase coexistences. 
It was shown there that the polydispersity can have a profound influence on the binodals 
of this system as compared with its monodisperse counterpart. One aspect which was 
not studied in ref.8 concerns the fact that, as a consequence of the Gibbs phase rule, a 
polydisperse system can also phase separate or "fractionate" into more than two fluid phases, 
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each phase differing both in average density and in size distribution. The study of such a 
polydispersity-induced equilibrium between three fluid phases is the main object of the 
present investigation. Meanwhile, we will also consider another question left unanswered 
in ref.8. This concerns the specific form given to the parent-phase size-distribution. In 
ref.8 we used some well-known expressions (Schulz-Zimm and log-normal) for this input 
size-distribution. These theoretical distributions describe particles with sizes ranging from 
zero to infinity. Since both the very small and very large particles are absent from the 
experimentally studied systems, one may wonder whether their presence in these theoretical 
expressions could not lead to artifacts in the phase diagram. The influence of these tails (for 
both small and large sizes) of the theoretical distributions on the coexistence properties of 
ref.8 will hence be investigated here before proceeding to the three-phase equilibria. 

In section [TI| we recall the vdW-description of polydisperse fluids as used in ref.8. Sec- 
tion H| is devoted to the influence of the tails of the parent-phase distribution on the results 
of ref.8. In section |IV] we show that the present vdW-description allows for the phase equi- 
librium between three fluid phases and we study the relative stability of the two-phase and 
three-phase equilibria in section |V|. Our conclusions are given in the final section [VT]. 



II. THE POLYDISPERSE VDW FREE-ENERGY AND ITS CRITICAL POINTS 



The thermodynamic properties of a fluid of spherical particles with a size-distribution 
can be determined from the following generalization of the vdW free-energy (see ref.8 for 
details): 

f(T, [p]) = k B Tj dap(a){H A3{ ^ a) ) - 1} 

+ ll d(T l d °'V((T,o')p((T)p(o') (1) 

where, f(T, [p]), is the free-energy per unit volume at the temperature T (kg being Boltz- 
mann's constant) for a fluid for which the number density of species a is p(c), while A(cx) 
denotes the thermal de Broglie wavelength of species a. Any functional dependence on the 
density distribution p(a) is, as usual, indicated by [p\. In (|l|), the excluded volume repulsions 
are represented by: 

E\p] = l- Jdav(a)p(a) (2) 



where, v (a) = ^ (_R(<t)) 3 , is the volume of a spherical particle of radius R(a), whereas the 
attractions are represented by: 

V(a,a') = J drV A (r;a,a'). (3) 

with Va{t] <j, cr') being the potential of attraction between two particles of species a and cr', a 
distance r = |r| apart. In the above a represents both a species label and the dimensionless 
polydispersity variable, R(a)/R(l), R{1) being the radius of the reference species, o = 1. 
The pressure, p(T, [p]), corresponding to (|l|) is: 
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p(T, [p]) = / dap(a) + - J da J da' V(a, a')p(a)p(a'). (4) 
whereas the chemical potential of species a, p(a,T, [p]), reads: 

+ J da'V(a,a')p(a'). (5) 

In the above, all integrals over <r extend over the whole domain for which p(a) is non-zero. 
As in ref.8, we write, p(a) = ph(a), where p = J dap(a) is the average density and h(a) 
the (normalized) size-distribution. The dimensionless average density will be written as, 
7] = pv(l), with v(l) the volume of the reference species a = 1, while the dimensionless 
temperature will be, t = fc B T/e(l,l), with e(l, 1) = — V(l, l)/8v(l) and V(l, 1) being the 
integrated amplitude of the attraction between two reference particles (cf.(|3|) and ref.8). 

In ref.8 we have considered several types of polydisperse interactions characterized by 
two parameters {l,n}, viz. v(a) = v(l)a 3n and V(a, a') = [(a n + a' n )/2f(aa') l V(l, 1). 
As shown there the model with {I — l,n = 0} has a phase behavior which is similar 
to that with the full polydispersity {/ = l,n = 1} but is simpler to study because its 
excess free-energy involves fewer moments, = / daa k h(a), of the size-distribution h(a). 
Henceforth, we will therefore consider only this {I — l,n — 0}-model, so as to simplify the 
calculations. In ref.8 we also considered several types of parent-phase size-distributions h (a) 
but found no qualitative differences between them. In order to avoid too much duplication 
of results, we will henceforth consider only the log-normal parent-phase size-distribution, 
ho(a) = c exp[— a ln 2 (a/b)], where the three parameters {a, b, c] are determined by imposing 
the first three moments of ho(a), = daa k ho(a), to be such that: 1) h^{a) be 

normalized (m^ = 1), 2) the average value of a be one, i.e. equal to the reference species 
(mj '* = 1), 3) the variance of h (a) be - = I — 1 > 0, with / being the polydispersity 
index (m^ = I)- In terms of / this yields: a = 1/2/n/, b = J -3 / 2 , c = I /^/2-KlnI and 

m (0) _ jfc(fc-l)/2_ 

When / = 1 (or a — oo) the system is monodisperse (h (a) = S(a — 1)) and its phase 
diagram consists of the usual vdW-binodal ending in the vdW critical point (r) c — 1/3, t c — 
32/27). For a modest polydispersity, say / = 1.02 or a = 50, the phase diagram is similar 
to that of Fig.4 of ref.8. For each parent-phase density rjo, there now is a different binodal. 
Each binodal is truncated upwards at a maximum temperature, t m , with the corresponding 
densities, r]i(t m ) and ry 2 (t m ), lying respectively on the cloud-point and shadow curves. For 
a critical value of r]o, t]q = r] c , the corresponding binodal passes through the intersection of 
the latter two curves. This occurs for t = t c and, since, rji(t c ) = 772 (i c ), the point (r] c ,t c ) 
is a critical point where the two coexisting phases (1 and 2) become identical. However, 
in contradistinction with the monodisperse case (1 = 1) the critical temperature is not the 
borderline between the one-phase and two-phase regions. Indeed, in the polydisperse case, 
the two-phase region extends to t' m > t c , where t' m corresponds to the maximum of the cloud- 
point and shadow curves. At the same time the critical point loses one of its attributes, i.e. 
for / 7^ 1, there also appears a second high-density critical point. The latter is polydispersity- 
inducedi but, as will be shown below, not necessarily thermodynamically stable. When the 
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polydispersity index / is increased still further the ordinary low-density (vdW) critical point 
moves to a higher density and a lower temperature, while the polydispersity-induced high- 
density critical point moves to a lower density and higher temperature, until for a limiting 
value of /, say 1 = 1*, the two critical points merge as illustrated in Fig.l. For I > I* there 
are no critical points and the phase behavior loses all contact with its monodisperse (1 = 1) 
counterpart. 

III. TRUNCATED PARENT-PHASE SIZE-DISTRIBUTION 

For a polydisperse system the input data for the study of its phase behavior involve, 
besides the temperature (t) and the average density (rj ), the parent-phase size-distribution 
h (a) (as in ref.8 the subscript zero refers to the parent-phase). In theoretical work the values 
of a are usually distributed over the full interval of all possible sizes, < a < oo, whereas in 
the samples used in the laboratory these sizes usually span a continuous but finite interval, 
a* < o < a**, with < a* < 1 and 1 < cr** < oo if o = 1 represents the reference size 
(usually taken to be the average size). As stated already in ref.8, the small-size (0 < o < a*) 
and the large size (a** < o < oo) tails of ho(a) will be harmless for strongly peaked 
distributions, i.e. for small polydispersities, but, as stated in the Introduction, they could 
lead to artifacts for large polydispersities characterized by broad size-distributions. Before 
studying the phase behavior of our system in the region of the (large) limiting polydispersity, 
/ I*, where any contact with the monodisperse system is lost, we will first assess the 
influence of these tails on the results of ref.8. To this end we again consider the log-normal 
distribution of section [TT] but normalize it now over the finite interval (a*, a**) and rescale 
a in such a manner that the average value of o computed over this finite interval be again 
equal to one, i.e. the reference particle remains unchanged. In Fig. 2 we show the effect 
of this truncation of the cx-domain on the form of ho(a) for some representative values 
of {a*, a**}. In Table 1 we show its influence on the thermodynamic data. The overall 
effect of this truncation is hence to render the system less polydisperse. It also shifts the 
coexistence densities to lower values. Since some of our coexisting phases have densities for 
which the fluid phases could become metastable with respect to solid phases, a truncation 
of ho(a) could render these fluid phases their thermodynamic stability. In any case, a 
quantitative comparison with experimental results will require phase diagrams calculated 
for truncated size-distributions and hence require information about a* and a**. Since, 
however, no qualitative differences (or artifacts) are found, we will continue henceforth with 
the "untruncated" h (a) log-normal distribution, as defined at the end of section || . 

IV. THREE-PHASE EQUILIBRIA: LOCAL STABILITY 

As seen in section [TI], for a sufficiently large polydispersity index (/ > /*) the phase 
behavior of our system loses any contact with its monodisperse counterpart. This is hence a 
good region to look for novel features such as three-phase equilibria. In the present section 
we will investigate how an initial (untruncated) log-normal parent-phase size-distribution, 
h (a), fractionates into (two or three) daughter phases when its polydispersity index / lies 
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on either side of the limiting value I* (I* = 1.072 or a* = 13.72) where the critical points 
disappear. We will focus on two cases: a = 15 or / = 1.066 and a = 13.5 or I = 1.074. 



A. Two-phase equilibria 

As discussed in detail in ref.8, these can be obtained by solving the two-phase coexistence 
conditions: 

M T > M) = p( t > N); li{<r, T, [pi]) = /i(a, T, \p 2 \) (6) 

where {pi(cr) = p\h\(a), p 2 (c) = P2h 2 (cr)} are the density distributions of the two daughter- 
phases resulting from the fractionation of the parent-phase, po(&) = Po^-o( cr )- The integral 
equation resulting from (^) can be solved as explained in ref.8. For the large /-values 
considered here we find two sets of binodals, each one being linked to a different critical 
point (see section |TJ). Since at these critical points the system is only marginally stable, i.e. 
we have S 2 f(T, [p]) = 0, we first inquire for the local stability (5 2 f > 0) of these solutions 
with respect to an infinitesimal change Sp(a) of p(a) = p n (a) {n = 1,2): 

^ = /^/^|^) M<Tl)M<T2) >0 (7) 

(see ref.8 for the explicit form of S 2 f). It turns out that, for the r] - values investigated, one 
of the two sets is locally unstable, i.e. it does not satisfy (0) (see Fig.3). 



B. Three-phase equilibria 

Similarly, the phase equilibria between three fluid-phases are governed by the solutions 

of: 

p(T, [ Pl ])=p(T,[p 2 ])=p(T, [p 3 ]) 

K a > T > M) = M ff > T > N) = K a > T > N) (8) 

where {pi(cr) = p\hi(a), /^(c) = P2^2(o"), Pz{.&) = Psh^a)} are the density distributions 
of the three daughter-phases. To solve (^) we can generalize the method used in ref.8. 
The h n (cr) (n = 1,2,3) are constrained by the particle number conservation of species 
a, X\hi{a) + x 2 /i2(cr) + ^3^3(0") = h (a), whereas the total particle number conservation 
implies X\ + x 2 + £3 = 1, x n being the fraction of the parent-phase particles which went 
into phase n (see ref.8 for details). The p n are constrained by the conservation of the total 
volume, X1V1 + x 2 v 2 + x 3 v 3 = v , with v n = l/p n {n = 0,1,2,3). This leaves us with 
two independent distributions, say h\{a) and h 2 (a), two independent number fractions, 
say xi and x 2 , and two independent number densities, say p\ and p 2 . Given [x\^x 2 ) and 
(pi,P2), the two distributions, h\{a) and h 2 (a), can be found by solving the system of two 
integral equations resulting from p(a,T, [pi]) = p(cr,T, [p 2 ]) and p(a,T, [p 2 ]) = p(a,T, [^3]). 
As in ref.8, these two equations can be rewritten as: h\{a) = ho(a)Hi(a) and h 2 (a) = 
ho(a)H 2 (cr), so that the normalization of the h n (a) can be expressed as, 1 = / dah {a)Hi{a) 
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and 1 = / dah (a)H2(cr), whereas solving the latter two equations together with, p(T, [pi]) = 
p(T, [p 2 ]) and p(T, [p 2 ]) = p(T, [p 3 ]) will determine (xi,x 2 ) and (p!,p 2 ). The full solution 
of a three-phase equilibrium problem for a polydisperse fluid is thus a rather complex task 
(hereby justifying some of the simplifying assumptions introduced above). In the case where 
the excess free-energy depends only on a finite number of moments of the distribution h(a), 
as is the case for the vdW free-energy, the above procedure can be simplified as the integral 
equations can then be transformed into a finite set of non-linear relations between these 
moments, as explained in ref.8. We found that the amount of labor involved is similar to 
that of the projection method of ref.6 but contrary to the latter, the present method involves 
no approximation to the basic equations (|^) or ([8]). In order to speed up the convergence 
of the solution method one may use the Powell-algorithmlll instead of the more traditional 
Newton-Raphson methodO used in ref.8. An example of a three-phase coexistence found in 
this way is shown in Fig. 4. We have verified that the case shown is locally stable, i.e. the 
p„0) (n = 1,2,3) satisfy (0). 

V. RELATIVE STABILITY OF THE TWO-PHASE AND THREE-PHASE 

EQUILIBRIA 

As seen in the previous section, in strongly polydisperse systems it is possible to obtain, 
for the same input data, locally stable two-phase as well as three-phase solutions. This then 
raises the question of how to separate the stable from the metastable transitions, i.e. of 
the global stability of these solutions. All the globally stable states belong to the convex 
envelope of the free-energy surface. This envelope is such that the tangent plane through 
any of its points never cuts the free energy surface. Of course, in the polydisperse case these 
"surfaces" are defined in the infinite dimensional functional space supporting the density 
distributions, p(o"). As a consequence, if p e (cx) belongs to the convex envelope of f(T, [p]) it 
must satisfy: 



for any p(cr). It is of course not possible to verify this infinite number of conditions but, in 
practice, it will suffice to verify (||) for those states p e (c) which are candidates for a stable 
equilibrium, i.e. for the solutions of eqs. (|6]) or (||), and for those states p(a) for which the 
distance between the free-energy surface and the tangent plane through p e (a) is extremal, 
i.e. for the p(cr) solution of: 



sible that there exist (in the same parameter region) higher-order equilibria involving four, 
five, etc. phases which could still invalidate @, but in view of the tremendous complexity 
of the global stability problem, we will limit ourselves here to the relative stability of the 
two-phase and three-phase equilibria found thus far. As a result of this search we find that 
the locally stable two-phase equilibrium solution of (H) is globally stable between an upper 






course still pos- 
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(t + ) and a lower (t_) temperature. Above the upper temperature (t > t + ) the parent-phase 
is stable, while for t < t + the parent -phase first fractionates into two phases while for t < £_ 
the system further fractionates into three phases. Two examples of such phase diagrams are 
shown in Fig. 5 and Fig.6. 



VI. CONCLUSIONS 

The investigation of the phase behavior of the polydisperse generalization of the vdW 
free-energy which was started in ref.8 for relatively modest polydispersities (/ ~ 1.02) has 
been extended here to larger polydispersities (I ~ 1.07). It is seen that when the polydisper- 
sity is increased the system's phase behavior gradually loses contact with its monodisperse 
counterpart. Indeed, for small polydispersities, the system investigated has two critical 
points, one of which is the polydisperse generalization of the vdW critical point, the other 
being polydispersity-induced. When the polydispersity is increased these two critical points 
merge for a limiting polydispersity (I* ~ 1.072). Above this threshold value (J*) there are 
no critical points. Similarly, when the polydispersity is increased, the two-phase region splits 
into two two-phase regions for the higher temperatures whereas for the lower temperatures 
a polydispersity-induced three-phase region appears. In order to determine the boundaries 
between these different regions of the phase diagram it is essential to investigate both the 
local and global stability of a large number of possible phase transitions. While some of these 
transitions are not even locally stable others are locally and globally stable. The situation 
quickly becomes very complex both physically and mathematically. 

The phase behavior investigated here should be of relevance to the study of colloidal 
dispersions of spherical particles for which the size-distribution is monomodal. A detailed 
comparison will of course require a further extension of theoresent investigation towards 
the solid phases, in particular their thermodynamic stability^ for the large polydispersities 
considered here. 
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Figure Captions 



FIG. 1. The low-density (circles) and high density (squares) critical points versus the 
polydispersity (I — 1 + ^). The limiting value of a below which the two critical points 
disappear is a* = 13.72 (I* = 1.073). Fig. la (resp. Fig. lb) displays r\ c (resp. t c ) versus a. 

FIG. 2. The full log-normal size distribution (dashed-dotted line) compared to a log-normal 
distribution truncated for a < cr* and a** < a (full line). Both distributions are normalized 
(m = 1 = itLq) and have the same average a- value (m x = 1 = m\). The case shown here 
corresponds to a — 15 (/ = 1.067) and a* = 0.72, a** = 1.33. 

FIG. 3. Two sets of binodals for a — 15 and i] = 0.68. The full-line binodals correspond 
to locally stable solutions of eq.(6) while the dotted-line binodals are locally unstable. Also 
represented (diamond) is the second (unstable) critical point. 

FIG. 4a A three-phase equilibrium for a = 15 and i] = 0.506 corresponding to a locally 
stable solution of eq.(8). The coexisting densities at t = 1.07 are r/i = 0.284, r] 2 = 0.652 and 
773 = 0.930. 

FIG. 4b Coexisting distributions hi(cr) (i = 1: full-line, i = 2: dotted-line and % = 3: 
dashed-dotted line) at t = 1.07 and for the a and f] values of Fig. 4a. 

FIG. 5 A phase diagram for a — 15 and i] = 0.506. The one-phase region corresponds to 
t > t + ~ 1.75, the two-phase region to t + < t < £_ ~ 1.24, while the three-phase region 
corresponds to t < 

FIG. 6 The same as Fig. 5 but for a = 13.5 and 770 = 0.506. Here t + ~ 1.75 and 

t- ~ 1.26. Note that while in Fig. 5 the new (third) phase is a low-density phase here it is 

an intermediate-density phase. 
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TABLES 





a** 


m 


V2 





oo 


0.24 


0.69 


0.30 


2.75 


0.24 


0.69 


0.48 


1.76 


0.23 


0.68 


0.56 


1.54 


0.20 


0.66 


0.70 


1.32 


0.15 


0.63 



TABLE I. Typical shifts in the coexisting densities of the low-density (rji) and high density 
(772) phase for t = 1, a = 15 after truncation of the log-normal parent-phase size-distribution for 
a < a* and a** < a. The corresponding monodisperse results are rji = 0.10, 772 = 0.61. 
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Fig. la , Bellier-Castella et al., JCP 




Fig. lb , Bellier-Castella et al., JCP 



1.18 
1.14 




Fig. 2 , Bellier-Castella et al., JCP 




Fig. 3 , Bellier-Castella et al., JCP 




Fig. 4a , Bellier-Castella et al., JCP 




Fig. 4b , Bellier-Castella et al., JCP 




Fig. 5 , Bellier-Castella et al., JCP 




Fig. 6 , Bellier-Castella et al., JCP 




